Inhomogeneous Ferromagnetism and Unconventional Charge Dynamics 
in Disordered Double Exchange Magnets 
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We solve the double exchange model in the presence of arbitrary substitutional disorder by using a 
self consistently generated effective Hamiltonian for the spin degrees of freedom. The magnetic prop- 
erties are studied through classical Monte Carlo while the effective exchange, Dij, are calculated by 
solving the disordered fermion problem, and renormalised self-consistently with increasing temper- 
ature. We present exact results on the conductivity, magnetoresistance, optical response and 'real 
space' structure of the inhomogeneous ferromagnetic state, and compare our results with charge dy- 
namics in disordered Lai-^Sr^MnOa. The large sizes, O(10 3 ), accessible within our method allows 
a complete, controlled calculation on the disordered strongly interacting problem. 
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The double exchange (DE) model provides the simplest 
example of strong coupling itinerant electron ferromag- 
netism, and is at the heart of the remarkable magneto- 
transport [1] in the manganites. The qualitative mag- 
netic features of the model are reasonably well under- 
stood now [2], although transport properties have been 
much less explored. The problem is more difficult in the 
presence of quenched disorder, which arises inevitably in 
real materials [3-7] from ionic substitution and resulting 
bond distortions. Disorder leads to an inhomogeneous 
magnetic state and enhances the density of low energy 
spin fluctuations. The resistivity arising out of structural 
disorder and magnetic scattering is usually large, compa- 
rable to the Mott resistivity [5,6], violates Mathiessens 
rule, and cannot be accessed by standard transport the- 
ory. The optical response reveals a strongly non Drude 
character [6,7], implying unconventional charge dynam- 
ics, and rapid loss in low energy spectral weight with 
rising temperature. 

Most of these features, which depend explicitly on 
the inhomogeneous spatial character of the magnetic and 
electronic state cannot be captured within 'mean field' 
approximations, including 'dynamical mean field the- 
ory' (DMFT), and current 'real space' approaches are 
severely size limited in three dimension. In this paper 
we use a new [8] Monte Carlo (MC) technique, which 
handles the interplay of disorder and spin correlation es- 
sentially exactly, to provide the first controlled results 
on charge dynamics in the disordered double exchange 
(DDE) model. We map out the phase diagram, clar- 
ify the nature of the inhomogeneous ferromagnetic state, 
and provide results on magnetotransport and optical re- 
sponse in the model. Our results are directly relevant to 
the 'coherent to incoherent' crossover [5-7] in disordered 
Rl-aSr^Mni-zAlzOa, and the properties of disordered 
metallic ferromagnets [9] in general. This approach can 
be readily extended to include phonon degrees of free- 
dom, of key importance in manganite physics. 



The DDE model has been proposed earlier to explain 
the ferro-metal (FM) to para-insulator (PI) transition, 
and the associated colossal magnetoresistance (CMR) in 
manganites. It is now accepted that electron-phonon cou- 
pling, and possibly phase coexistence, are important [1] 
in understanding CMR. While the DDE model can ex- 
hibit a FM-PI transition, at large disorder, in real ma- 
terials even weak electron-phonon coupling will strongly 
affect the 'Anderson localised' PI phase. We will illus- 
trate the FM-PI transition, but focus more on the 'metal- 
lic' phases, of relevance to the disordered Rei-^Sr^MnOs 
family at x > 0.3. In these materials phonon effects are 
not significant and the effect of spin fluctuations on dif- 
fusive charge dynamics can be directly probed. 

We study the following model: 
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-t are nearest neighbour hopping on a cubic 



lattice and the on site disorder is distributed uniformly 
between ±A/2. We set Ju/t — > oo. The parameters in 
the problem are A/t, and density n (or chemical potential 
fi). We set t = 1, fixing our basic energy scale, and 
assume the core spins to be classical, with |Sj| = 1. 

There are two key steps in solving for the magnetic 
and transport properties of a model like this: (i) Eval- 
uate the correlated spin distribution, P{Si}, controlling 
the magnetic response, by 'integrating out' the electronic 
degrees of freedom. {Si} denotes the full spin configura- 
tion, (ii) Solve for charge dynamics; resistivity, optical 
response, etc, in the background of structural disorder 
and equilibrium spin configurations. 

Previous studies of the model have used variational 
mean field (VMF) theory [10] and DMFT [11,12] to ac- 
cess the magnetism. Transport properties have been 
analysed within DMFT [11] and by using 'scaling the- 
ory' [13] in the limits of a spin polarised (T = 0) state 
and a fully spin disordered (T ^> T c ) state. Finally, real 
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FIG. 1. (a). Approximate 'universal' phase diagram of the 
DDE model. For a specified electron density, n, A is measured 
in units of the T = critical disorder, Af M (n), while T is 
measured in units of T c (n, A = 0). Our data for n = 0.1 
(circles) and n = 0.3 (squares) are described approximately 
by a common fit. (b). Actual phase diagram for n = 0.3 The 
'boundary' separating the PM and PI phases is a crossover. 

space MC technique has been used [14] to study the mag- 
netic transition and the spin wave spectrum. Both VMF 
and DMFT 'factor' the correlated spin distribution into 
a product of single site distributions. This loses out on 
spin correlations and trivialises the paramagnetic phase. 
In addition, DMFT misses out key vertex corrections in 
the 'disorder average' process, losing electron localisation 
effects. The 'scaling theory' approach to electronic trans- 
port is exact for T — ► and T 3> T c , but not of much 
use in the crucial regime around T c . The only approach 
which implements step (i) exactly is Monte Carlo, but 
the computational cost of these simulations limits avail- 
able sizes to ~ 4 3 in three dimension [15]. The small 
linear dimensions, and the large finite size gaps, make 
it impossible to reliably estimate transport properties. 
This is where our method, below, allows a breakthrough. 

For Jn/t — ► oo, a standard transformation and projec- 
tion [8] can be used to map on H to a 'spinless fermion' 
problem with core spin dependent hopping amplitudes: 
H = -*£<«> fij( e^nh +h.c) + £.( Ci - rim. The 
hopping amplitude g^ — fije l ^ ij , between locally aligned 
states, can be written in terms of the polar angle 



and azimuthal angle 



>i) of the spin Si as, cos^-cos-f 
The 'magnitude' of the over- 



lap, fij = + Sj.Sj)/2, and the phase is specified by 
tan$ij = lm{g ij )lRe{g ij ). 

The fermions in this 'quadratic' problem move in the 
background of quenched disorder ej and 'annealed disor- 
der' in the hopping amplitudes gij. To exploit the 'non 
interacting' character of the fermion problem we need 
to know the relevant {/, 3>} configurations, controlled by 
H e ff{f,&} = -^logTre^PH . The corresponding Boltz- 

mann distribution is P{f, $} oc e~ Heff ^'^. Our key 



model, is determined self consistently as the thermal 
average of the 'mixed' spin-fcrmion operator, = 
(e** w 7j-7j + h.c) over the assumed equilibrium distribu- 
tion. The self-consistent Dy are solved for via MC for 
a specified /i, {e^}, and T. The can be spatially 
strongly inhomogeneous and also significantly tempera- 
ture dependent. At self-consistency fermionic averages 
are computed over equilibrium spin configurations. We 
work at constant n, fixing fi through iteration. Since 
we avoid the expensive spin update procedure of 'exact' 
MC and use diagonalisation only to compute the we 
can access sizes ~ 10 3 , compared to 4 3 in the standard 
approach. Transport properties, at equilibrium, are com- 
puted exactly using the Kubo formula, employing sizes 
~ 8 3 — 10 3 . We systematically check for size dependence 
in our transport results. The conductivity results are in 
units of (ire 2 )/hao, the Mott limit being (0.03e 2 )/?iao. 

Fig.l shows the "global" phase diagram of the DDE 
model. We have studied the problem at n = 0.3 and 
n = 0.1, varying disorder from the perturbative end to 
the localisation regime. Fig.l. (a), superposes the results 
at n — 0.3 and n = 0.1, appropriately scaling the disor- 
der and temperature (see caption). The critical disorder 
for the PM-PI crossover, or the disorder dependence of 
T c need not be 'universal' but seems to follow the same 
overall trend at moderate n. Fig. 1.(6) shows the 'true' 
phase diagram specifically at n = 0.3. The T c of DE 
models is approximately related to the internal energy 
change: T c log(2S + 1) - £{T C ) - £({)), which in turn is 
related to the kinetic energy, K, at T = 0. At small A 
and T = 0, K(A) - K (0) - x d A 2 , where \d is the 'lo- 
cal' density response function. At strong disorder, in the 
'localised' phase, K(A) oc i 2 /A. The correspondence of 
these limits with the inferred T C (A) is visible in Fig. 1.(6). 
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FIG. 2. (a). Temperature dependence of magnetisation at 
n = 0.3, with A increasing from — 24. The inferred T C (A) is 
shown in Fig.l. (6). Comparison of m(T) obtained via exact 
MC and -He// on a 4 x 4 x 4 system at n = 0.3. Disorder 
average over 8 realisations in both case. (c). The scaled mag- 
netisation m(T/T c (A)), at n = 0.3,for A = and A = 24. 
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FIG. 3. Real space images, top surface of a 10 x 10 x 10 cube 
with PBC, at n = 0.3 and A = 6, for a specific realisation 
of disorder. First row: electron density distribution, (n(r))T, 
second row: effective exchange, D(r, r + 8 : T), third row: 
nearest neighbour spin correlation, f 2 = (S(r).S(r + S)) T - 
Temperatures along the row are, T/T c ~ 1.2, 0.9, 0.6 and 0.2. 
Dark regions correspond to high density (top row), strong 
exchange (central row), and strong FM correlation (bottom). 

Note that even at strong localisation, A <~ 24, local 
excursion of the electrons can still sustain a T c ~ 0.05, 
which would be - 75K for t ~ 150meV. DMFT and VMF 
also provide qualitatively similar trends in T c . 

Fig. 2. (a) shows m(T) with increasing disorder, at n = 
0.3, from which the data in Fig. 1.(6) was inferred, while 
Fig. 2. (6) compares the m(T) computed with H e ff with 
the 'exact' MC at A = 6. Although the asymptotic be- 
haviour of T C (A) is easy to motivate, Fig.2.(c) provides 
the first indication that the properties of the disordered 
ferromagnet cannot be understood by merely scaling T c . 
Even though the T = state is fully polarised, at low 
finite T the 'weak' bonds lead quickly to local disorder- 
ing. The magnetisation m{T/T c ) falls more sharply at 
low temperature in the disordered system [9] . The inho- 
mogeneous character, and correspondingly wide distri- 
bution, of the 'exchange' leads to an overall increase in 
the density of low energy magnetic excitations. In addi- 
tion to suppressing m(T) it enhances the specific heat, 
Cv(T/T c ), and reduces the spin wave stiffness Dsw- 

The real space images, Fig. 3, illustrate the inhomoge- 
neous freezing into a ferromagnetic state in a 'cooling se- 
quence', left to right. The density profile (n(r)), first row, 
remains unchanged over the T range of interest, 0— 1.2T C . 
The exchange, Dij(T), central row, is strongly spatially 
inhomogeneous, but as a whole not strongly T depen- 
dent. The bond distribution, P(D,T), reveals [16] that 
weak bonds at low T get quickly weakened with increas- 
ing T while strong bonds are essentially T independent. 
The nearest neighbour spin correlation, /2, highlights the 
inhomogeneous disordering of the polarised state with in- 
creasing T (white regions, two panels on the right), and 
the surviving local order for T > T c (left panel). 

These results suggest the possible correlation between 



spatial inhomogeneity in the magnetic state and the bulk 
thermodynamics. The interplay of thermal spin disorder 
and 'frozen' structural disorder also affects the charge 
dynamics. Increasing A increases the residual resis- 
tivity, with po <~ PMott at A ~ 10. The resistivity 
in the paramagnetic phase (p» n /), tracked at T = 0.4 
in Fig. 4. (6), is not simply the additive contribution of 
structural and spin disorder. If Mathiesscns rule were 
obeyed, pi n j — po should have been constant. The de- 
viation arises from interference between structural and 
magnetic scattering and is clearly observed in the metal- 
lic Rei_ x Sr x Mn03, at x = 0.4, with Re being, La, Pr, 
Nd [6], and in Lai_ a; Sr 2 ;Mni_2Al z 03, at x — 0.3, vary- 
ing z [5]. This 'interference' is beyond Boltzmann the- 
ory. With a <~ 4A, pMott ~ 5 mCicra in the manganites. 
For the 'disordered' LaSr family, as po/ PMott varies from 
0.005 — 0.04 with increasing disorder, the corresponding 
pinf I PMott increases from 0.4 — 1.5 [5,6]. 

Fig. 4. (a) shows the normalised resistivity Sp(T) = 
(p(T) - p(0))/(p(0.4) - p(0)) for A = 0,2,6. The shift 
in T c in this disorder regime is quite small, but Sp(T/T c ) 
quickly changes character with increasing disorder. The 
more prominent short range spin fluctuations in the dis- 
ordered system couple to the diffusive electrons leading 
to a sharper rise in p(T). While most of the rise in the 
'clean' system occurs in the vicinity of T c , the rise is 
spread over a wide interval in the disordered system. 
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4. (a). Temperature dependence of resistivity at 
n = 0.3. We plot p(T,A) - p(0, A) normalised by 
p(0.4, A) — p(0, A), shifting successive curves by 0.5 for clar- 
ity. T c marked by an arrow on each curve, (b). Resistivity 
at T = 0, po, and at T = 0.4, Pi„f, with increasing dis- 
order. Note the clear violation of Mathiessens rule even at 
A = 2 (see text), (c) Resistivity normalised to maximum 
value (~ 10 4 ), and m(T) at n — 0.1 and A = 12 illustrating 
the FM — * PI transition, (d). Variation in T c with residual 
resistivity. Note the logarithmic x scale, (e). The magnetore- 
sistance (p(T,0) - p(T, h))/p(T, 0) at n = 0.3: A = (lower 
set) and A = 6 (upper set). The A = 6 set has been vertically 
shifted by 0.8. Field values are h = 0, 0.02, 0.05, 0.10. 
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FIG. 5. Optical conductivity at n — 0.3. (a). A = 0, 
(6). A = 4, (c). A = 6, at h = 0, and (d). A = 6 at 
T ~ T c . Other parameters indicated in the figure. Inset in 
panel (d) shows the integrated spectral weight, N e ff(T), (see 
text) derived from the data in panel (c). Data averaged over 
4 — 8 realisations of disorder. 



Fig.4.(c). shows data at n = 0.1 and A = 12 to illus- 
trate the FM-PI transition. We work with this lower den- 
sity because the MIT is easier to access. We solve the full 
problem for sizes 6 x 6 x L, with L = {6, 12, 24}, disorder 
average, and extrapolate the computed 'd.c conductivity' 
to L — > oo. This is crucial to capture 'Anderson locali- 
sation' in finite systems. The transport in the insulating 
phase is controlled by activation to the mobility edge. 

The correlation between T c and residual resistivity has 
been experimentally explored [3,4] and Fig.4.(e) high- 
lights the fall in T c as po heads towards the Mott limit. 

Our results on MR, Fig.4.(e), indicate that in the in- 
termediate disorder regime the MR is actually smaller in 
the disordered problem than in the clean system. The 
MR is 1 - p(T, h)/p(T, 0). The minimum p(T, h) is lim- 
ited by po, the residual resistivity. With increasing A, p 
grows, limiting the MR. The MR rises again only when 
A is large enough to drive a finite T, FM-PI transition. 

Finally, the optical response of the system, at n = 0.3 
and A = 0,4 and 6, is shown in Fig. 5. Panel (a), shows 
<j{oS) at A = and the Drude feature survives even for 
T <~ 2T C with Udc ~ 6<jmo«- For A = 4, however, panel 
(b). reveals that cr(u)) changes from a Drude form for 
T < T c to an essentially 'flat' incoherent response at 
T <~ 2T C . At 2T C , adc ~ ZcMott- This response is roughly 
like Lai-zSr^MnOs at x <~ 0.4 [6]. At even larger dis- 
order, A = 6, the response is non Drude even at T = 0, 
and becomes markedly so, with a finite w peak, as T 
is increased. The apparently t/uj rise results from the 
intimate coupling of the diffusive electrons to spin fluc- 
tuations via the hopping modulation. For A = 6, Ode ~ 
Zo-Mott at T = 0, and a dc ~ OMott at T — 2T C . This is 
like the response in (Pr,Nd)i_ a; Sr ;Z ;Mn03 at x <~ 0.4 [6]. 
We predict that the response will be similar at x = 0.4 



in Lai_ K Sr a; Mni_ z Al z 03 at z ~ 1% - 2%. The inset to 
panel (d). shows the integrated low frequency spectral 

weight N ef f(u',T) = J^' a(u,T)du, at oj' = 2. There is 
a 40% loss in spectral weight in the T range — 2T C . Our 
energy cutoff in N e ff roughly corresponds to w' ~ 0.3eV 
in the manganites. Panel (d). shows the magneto-optical 
response at A = 6 and T <~ 1.2T C . 

The non Drude relaxation with cr(u>) having a finite lo 
peak occurs in a regime where dp/dT > 0, as in conven- 
tional metals, but the charge dynamics is highly diffusive 
as expected in a system with strong 'effective disorder'. 
The change from 'coherent' to 'incoherent' dynamics oc- 
curs in a(u),T) when Pdc{T) ~ PMott- This general fea- 
ture is true of the PM phase of all the manganites. 

In conclusion, we have discussed the inhomogeneous 
magnetism in disordered double exchange magnets, and 
provided the first "exact" results on transport and optical 
response. The dependence of the non Drude relaxation 
on disorder and temperature is fully consistent with the 
metallic manganites. The extension of our method to 
include phonons will allow a complete solution of charge 
response in the manganites. 
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